Multi— Epoch VLBA Observations of the Compact Wind— Colhsion 
Region in the Quadruple System Cyg OB2 #5 



Sergio A. Dzib^, Luis F. Rodriguez^'^, Laurent Loinard^''^, Amy J. Mioduszewski^, Gisela 

N. Ortiz-Leon^ and Anabella T. Araudo^ 

s . dzibOcrya . unam . mx 

Received ; accepted 



^Centro de Radiostronomfa y Astroffsica, Universidad Nacional Autonoma de Mexico, 
Morelia 58089, Mexico 

^Astronomy Department, Faculty of Science, King Abdulaziz University, P.O. Box 80203, 
Jeddah 21589, Saudi Arabia 

■^Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

''National Radio Astronomy Observatory, 1003 Lopezville Road, Socorro, NM 87801, USA 



- 2 - 
ABSTRACT 

We present multi-epoch VLBA observations of the compact wind coUision 
region in the Cyg 0B2 #5 system. These observation confirm the arc-shaped 
morphology of the emission reported earlier. The total flux as a function of time 
is roughly constant when the source is "on" , but falls below the detection limit 
as the wind collision region approaches periastron in its orbit around the contact 
binary at the center of the system. In addition, at one of the "on" epochs, the 
flux drops to about a flfth of its average value. We suggest that this apparent 
variation could result from the inhomogeneity of the wind that hides part of 
the flux rather than from an intrinsic variation. We measured a trigonometrical 
parallax, for the most compact radio emission of 0.61 ± 0.22 mas, corresponding 
to a distance of 1.65 44 kpc, in agreement with recent trigonometrical parallaxes 
measured for objects in the Cygnus X complex. Using constraints on the total 
mass of the system and orbital parameters previously reported in the literature, 
we obtain two independent indirect measurements of the distance to the Cyg 
0B2 7^5 system, both consistent with 1.3-1.4 kpc. Finally, we suggest that the 
companion star responsible for the wind interaction, yet undetected, is of spectral 
type between BO. 5 to 08. 

Subject headings: stars: individual (Cyg 0B2 #5) — stars: winds — radiation 
mechanisms: nonthermal — astrometry 
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1. Introduction 

Non-thermal radio emission from a massive binary system is usually produced in the 
interaction of their strong winds (see De Becker 2007 for a recent review). An interesting 
case is the quadruple massive system Cyg 0B2 #5 (V729 Cyg, BD +40 4220) that is 
the only multiple system known so far to harbor two radio-imaged wind-collision regions 
(Ortiz-Leon et al. 2011). 

Cyg 0B2 ^5 is a radio-bright massive multiple system located in the Cygnus 0B2 
association, one of the most prominent massive young clusters known, containing over 100 
O stars. The most massive component of Cyg 0B2 ^^5 (i.e. the primary) is an eclipsing, 
contact binary consisting of two 0-type supergiants with a 6.6-day orbital period. We shall 
refer to this central object as Cyg 0B2 ^^5 A (containing the individual stars Cyg 0B2 ^^5 
Aa and Cyg 0B2 ^^5 Ab) . Contreras et al. (1997) suggested that a radio source located 
0"8 to the north-east of Cyg 0B2 ^^5 A is the shock formed by the interaction between 
the wind from Cyg 0B2 ^^5 A and that of a B-type star first reported by Herbig (1967) 
and located 0''9 to the north-east. We shall refer to this B star as Cyg 0B2 ^^5 C, and to 
the wind-collision region between this star and Cyg 0B2 ^^5 A as WCR(A-C). Recently, 
Ortiz-Leon et al. (2011) used the Very Long Base Array (VLBA) telescope to obtain a high 
angular resolution image (~10 mas) and reported the detection of a compact wind-collision 
region with possible non-thermal radio emission. This structure is formed by the interaction 
between the wind driven by Cyg 0B2 ^^5 A and that of an unseen companion (that we shall 
call Cyg 0B2 ^^5 B), an early B type star, with an orbital period of ~ 6.7-yr (Kennedy 
et al. 2010) Q. As this star orbits around Cyg 0B2 ^^5 A, the resulting WCR (hereafter 
WCR(A-B)) is undetected around what is presumed to be orbit periastron, yet can be 

^Kennedy et al. (2010) refer to Cyg 0B2 #5 Aa, Ab, B and C as Star A, B, C and D, 
respectively. 
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detected as the orbit progresses to apastron. Ortiz-Leon et al. (2011) noticed that the 
separation between the WCR(A-B) and the position of the contact binary is 12 mas, i.e. 
smaller than the nominal radius (23 ± 12 mas; Rodriguez et al. 2010) of the region where 
the optically-thick free-free radiation of the Cyg 0B2 #5 A binary system is produced. 
Therefore, the WCR(A-B) emission should in principle be undetectable due to the free-free 
opacity. To resolve this conundrum, these authors proposed that the wind driven by the 
Cyg 0B2 #5 is very inhomogeneous, as suggested for other early massive stars (Blomme et 
al. 2010; Muijres et al. 2011). However, the time variability could be due to other plasma 
effects (Pittard & Dougherty 2006; Pittard 2009). 

Being part of the Cyg 0B2 region, Cyg 0B2 #5 is expected to be located at the same 
distance. Traditionally, the Cyg 0B2 region has been beheved to be at a distance of 2.1 
kpc (Reddish et al. 1966), although it has, more recently, been argued to be somewhat 
nearer (at about 1.7 kpc; e.g. Massey & Thompson 1991). Even shorter distances have 
been proposed. For instance, Hanson (2003) suggested d ~ 1.4 kpc or even less. Direct 
trigonometric parallaxes were obtained recently by Rygl et al. (2012) who measured the 
distances to several methanol and water maser related to objects in the Cygnus X star 
forming complex. They found that most of them are consistent with a distance of 1.40 
± 0.08 kpc, and suggest that the Cyg 0B2 region is located at a similar distance. Using 
similar techniques, Zhang et al. (2012) measured the distance to the red hypergiant NML 
Cyg to be 1.61 ± 0.12 kpc and suggested that it hes on the far side of the Cyg 0B2 region. 
Thus, most recent, rehable measurements suggest a distance of 1.4 kpc for the Cyg 0B2 
region as a whole, and for Cyg 0B2 ^^5 in particular. Interestingly, however, Linder et 
al. (2009) obtained a distance estimate for Cyg 0B2 #5 itself of 925 ± 25 pc from a light 
curve analysis. A direct astrometric study to Cyg 0B2 #5 could help solve the discrepancy 
between the distance suggested by Linder et al. (2009) and the recently measured distances 
for other objects in the region. 
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In this paper, we present the analysis of a series of 12 VLBA observations of the 
WCR(A-B) in Cyg 0B2 ^^5, covering a total time span of 1.6 years. These data will be 
used to study the nature of the emission mechanism in the WCR(A-B), and the distance to 
Cyg 0B2 7^5. In the past, our team has successfully used multi-epoch VLBA observations 
to measure distances to low and intermediate mass stars with compact nonthermal radio 
emission (see Dzib et al. 2010 and 2011 for recent results). In the case of Cyg 0B2 ^^5 the 
wind collision region is resolved with VLBA observations, but should still be sufficiently 
compact to enable usable astrometry. A description of the observations, their calibration 
and imaging are described in Section 2. The results are presented in Section 3, where we 
also describe the structure and variability of the emission. Additionally a rough estimation 
of the trigonometric parallax is also presented in this section. The discussion of the results 
is presented in Section 4, and we finish with the conclusions in Section 5. 



2. Observation and data calibration 

Cyg 0B2 #5 was observed at A = 3.6 cm (z/ = 8.42 GHz) with the VLBA at 12 epochs 
between 2010 December and 2012 July. The first observation (in 2010 December) was 
designed as a detection experiment and was reported by Ortiz-Leon et al. (2011). Following 
this successful detection, we initiated a series of multi-epoch observations. The separation 
between successive observations in those multi-epoch data was about 1.5 months (Tabled]). 
They were usually realized at a recording rate of 512 Mbps but the epoch 2011 October was 
recorded using a rate of 2 Gbps, to test new VLBA equipment. As a consequence, the noise 
level for this epoch was significantly better. 

The observations consisted of series of cycles with two minutes spent on source, and 
one minute spent on the main phase-referencing quasar J0218+3851, located 3?6 away. 
Every 15 minutes, the radio bright X-ray binary Cyg X-3, located at 20' from Cyg 0B2 
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Fig. 1. — Images of Cyg 0B2 ^^5 obtained with the VLBA. The contours are at -3, 3, 4.5, 
6, 7.5, 9 etc. times the 3(T noise level in the image (see Tabled]). The synthesized beam is 
shown at the bottom left. 

^5, was also observed. Although Cyg X-3 is very close to Cyg 0B2 #5 it was not used as 
a primary calibrator due to its extreme variability in both flux density and morphology and 
the fact it is heavily scattered at the observing frequency (Desai & Fey 2001). As we shall 
see, however, it was used for secondary calibration. The total duration of the observations 
was 5 hours for the 2010 December observation and 2 hours for the other epochs. 

The data were edited and calibrated using the AIPS software, following standard 
procedures for phase-referenced VLBA observations. The calibration determined from 
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the observations of J0218+3851 was applied to both Cyg 0B2 #5 and Cyg X-3. At this 
point, most of the phase errors left are caused by the 3? 6 separation between the source 
and the phase calibrator. To remove most of this, Cyg X-3 was self-calibrated in phase and 
the incremental phase gains determined from that self-calibration were interpolated and 
applied to Cyg 0B2 #5. In most of the epochs, this resulted in a significant improvement 
in the quality of the image of Cyg 0B2 ^^5 compared to a direct phase transfer between 
J0218+3851 and Cyg 0B2 #5. For epochs 2011 August and 2011 December, however, 
this step worsened the image, presumably because of the high variability of Cyg X-3. For 
those epochs, we use the images obtained without applying this step. Cyg 0B2 ^^5 is quite 
resolved, so the best images were made by limiting the maximum uv length to 60000 kA 
and using always natural weighting. The total flux for each detection in these images was 
determined using the AIPS task IMSTAT. The rms noise level in the final images (shown 
in Figured]) was 43-110 /iJy beam^^ (Tabled]). Additionally, high resolution images were 
made using the whole uv range. The epochs at which the source was significantly detected 
in this second set of images are shown in Figure El 

As we will see, significant variations in the flux of Cyg OB27^5 were found. It is, 
therefore, worthwhile to examine in some detail the accuracy of the absolute flux calibration. 
We followed the standard VLBA recipes, based on gain and system temperature, for the 
absolute flux calibration. These methods are applied similarly to the target (Cyg 0B2 ^^5) 
and the gain calibrator (J0218-I-3851) so any systematic error should affect equally both 
sources. We inspected the obtained flux of J0218+3851 and found no systematic differences 
from epoch to epoch. Instead, we found a random scatter of about 0.035 Jy around a mean 
value of 0.308 Jy. Thus, the dispersion is about 10% of the mean and we will use that figure 
as our flux uncertaintyjj 



^Note that J2018+3851 is likely to be intrinsically moderately variable, and part of the 
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Given the 3?6 separation between the gain cahbrator and the target, part of the target 
flux might be lost under adverse weather conditions due to decorrelation. This can also be 
discarded using the observations of Cyg X-3 that were intertwined with the observations of 
Cyg 0B2 7^5. Cyg X-3 and Cyg 0B2 ^5 are very near each other (separation of about 
0°3) while the phase calibrator is 3°6 degrees away. Thus, if the phase transfer from the 
phase calibrator to the target were responsible for significant decorrelation on Cyg 062^^5, 
it should also induce significant decorrelation on Cyg X-3. We checked that this was not 
the case in any of our images, by comparing the flux measured in our VLBA observations 
with (publicly available) monitoring observations of Cyg X-3 made with the AMI-large 
array telescope at 15 GHzlfl 



3. Results 

3.1. Structure and properties of the emission 

From Figure [T] we can see that when the source is detected, the structure is similar to 
that reported by Ortiz-Leon et al. (2011): an arc-like morphology, with apex located to 
the west, extended mostly along the north-south direction. Following Ortiz-Leon et al. 
(2011), we attribute this morphology to the interaction between the strong wind driven by 
the contact binary Cyg 0B2 ^5 A located about 12 mas to the west of the WCR(A-B) (see 
Figure 2 in Ortiz-Leon et al. 2011), and the wind driven by an unseen companion (Cyg 0B2 

observed scatter might reflect this intrinsic variability. As a consequence, our accuracy might 
well be better than the quoted 10%. 

^ The AMI-large array telescope is supported by STFC and the Uni- 
versity of Cambridge and their monitoring can be accessed by the URL 
http : / / www. mr ao . cam. ac.uk/ telescopes / ami / index . html 
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Table 1: Observed epochs, flux densities (with statistic and systematic errors), and noise 
levels of the lo w resolution VLBA Images. 



Mean UT date 


Modified 




a 


(yyyy-mm.dd) 


Julian Day 


(mJy) 


(/xJy beam""^) 


2010.12.13 


55543.90 


3.00±0.13±0.30 


55 


2011.01.31 


55592.85 


0.66±0.10±0.07 


80 


2011.03.21 


55641.72 


2.98±0.18±0.30 


89 


2011.05.14 


55695.57 


3.51±0.16±0.35 


86 


2011.07.08 


55750.42 


3.24±0.16±0.32 


76 


2011.08.30 


55803.28 


2.57±0.13±0.26 


72 


2011.10.29 


55863.11 


1.66±0.07±0.17 


43 


2011.12.10 


55906.00 


2.58±0.13±0.26 


68 


2012.02.12 


55969.82 


1.33±0.16±0.13 


110 


2012.03.30 


56016.69 


0.98±0.14±0.10 


70 


2012.05.23 


56070.55 


<0.24 


80 


2012.07.31 


56139.36 


<0.25 


84 
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Fig. 2. — Images at higher resolution for the epochs when Cyg 0B2 ^^5 is still detected. 
The contours are at -3, 3, 4.5, 6, 7.5, 9 etc. times the 3cr noise level in the image (see Table 
[2]). The synthesized beam is shown at the bottom left. 

7^5 B) located just east of the apex (see again Figure 2 in Ortiz-Leon et al. 2011). This is 
also the preferred scenario of Kennedy et al. (2010) to explain the VLA observations of the 
system. In the high resolution images (Figure [2]), most of the flux is resolved out, and only 
the most compact radio emission remains detectable. The implied brightness temperature 
in these cases is in the range of 0.5-1.5 xlO^ K (Table [2]). 

The existing VLA observations of Cyg 0B2 ^^5 imply that the average flux from the 
WCR(A-B) during high state should be ~ 4.0 mJy at 8.4 GHz (as compared to the total 
flux of the stellar winds plus WCR(A-B) of ~9 mJy; Kennedy et al. 2010). Interestingly, 
however, the fluxes observed for WCR(A-B) in the VLBA observations reported here (Table 
H]) are significantly below the value of ~ 4.0 mJy (~ 3 mJy on average). This could be 
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taken to indicate that the flux in the current "on" state is systematically lower than the 
previous "on" states. To check this, we searched the archive of the Very Large Array (VLA) 
for recent observations of the system. We found an observation of the project 12B-165 
taken in 2010 August in C band and D configuration. The source was detected at 9.89 ± 
0.38 mJy, consistent with previous detection at the same band in the same configuration in 
previous "on" states (see Kennedy et al. 2010). We suggest, instead, that about 25% of the 
emission is resolved out by the VLBA. 

In Figure [3] the fluxes from Tabled] are plotted as a function of time. From December 
2010 to December 2011, the flux is roughly constant at about 3 mJy, except for the 
observation of 2011 January. After December 2011, the flux starts to decrease, and the 
source is finally not detected in 2012 May and 2012 July. This flux behavior is consistent 
with that expected from previous observations of Cyg 0B2 ^^5, and qualitatively in 
agreement with the best formal model (i.e. s=0) of Kennedy et al. (2010). Quantitatively, 
however, the best model of Kennedy et al. (2010) predicted the flux to drop about 7.5 
months later than it actually did (see Figure [3]). 

3.2. Astrometry 

For astrometric studies, point-like objects are preferred, and the WCR(A-B) in Cyg 
0B2 7^5 is clearly not the ideal type of target. However, the six high- resolution images 
shown in Figure [2] can be used to perform some astrometry, albeit with less accuracy than 
would be achievable for very compact sources. In particular, the extension of the source in 
declination induces large positional errors along that direction. Although it does include 
two successive epochs of maximum parallactic elongation (March and September), the total 
time span covered by our observations is only 8.8 months, and is also not optimal. Bearing 
in mind these limitations, we attempted an astrometric study of the WCR(A-B) in Cyg 
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Fig. 3.— Time evolution of the total flux of the WCR(A-B) in Cyg OB2#5, from the 
low resolution VLBA images. The upper limits are at Scr. The dotted curve shows the 
expectation from the best fit of Kennedy et al. (2010; the s=0 model) scaled by 0.75 to 
account for the missing flux (see text), while the sohd curve shows the same fit offset by 7.5 
months. 

0B2 #5. 



Table 2: Calendar dates and their corresponding Julian days, measured source positions, flux 
densities, noise levels, and brightness temperatures of the high resolution VLBA images. 



Mean UT date 


Julian Day 


a(J2000.0) 




<5(J2000.0) 




U 


(T 


Tb 


(yyyy.mm.dd/hh;mm) 




20''32™ 




+41°18' 




(mjy) 


(/ijy bcam^^) 


(lO'^ K) 


2011.03.21/17:17 


2455642.22 


22?421114 


OfOOOOlS 


18'.' 89983 


0'.' 00042 


2.17±0.43 


82 


1.24 


2011.05.14/13:44 


2455696.07 


22? 421094 


0?000012 


18'.' 90079 


0'.' 00039 


2.57±0.42 


75 


1.55 


2011.07.08/10:08 


2455750.92 


22? 42 1036 


0?000013 


18'.' 89937 


0'.' 00034 


1.32±0.27 


76 


0.56 


2011.08.30/06:40 


2455803.78 


22? 420999 


0? 000018 


18'.' 89792 


0'.' 00046 


1.37±0.32 


81 


1.02 


2011.10.29/02:44 


2455863.61 


22? 420940 


0?000012 


18'.' 89689 


0'.' 00042 


1.40±0.22 


41 


0.74 


2011.12.10/23:55 


2455906.50 


22? 420949 


0?000015 


18'.' 89466 


0'.' 00050 


2.19±0.44 


68 


0.89 
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The positions of the source in our VLBA data were determined from the high resolution 
images shown in Figure [2] using a two-dimensional Gaussian fitting procedure (task JMFIT 
in AIPS) and are given in Table [2J JMFIT provides an estimate of the position error based 
on the expected theoretical astrometric precision of an interferometer (Condon 1997); these 
errors are quoted in Columns 4 and 6 of Table El To obtain the astrometric parameters from 
these data, we used the single value decomposition fitting scheme described by Loinard et 
al. (2007). The necessary barycentric coordinates of the Earth, as well as the Julian date 
of each observation, were calculated using the Multi-year Interactive Computer Almanac 
(MICA) distributed as a CD ROM by the US Naval Observatory. The reference epoch was 
taken at the mean of the four observations: JD 2455774 = 2011.6. The best fit is shown in 
Figure HJ and yields the following astrometric elements: 

aj2oii.6 = 20^32"22?421021 ± 0.000006 

5j2mi.& = 41°18'l8"89823± 0.00032 

fj^a cos 5 = —1.64 ± 0.98 mas yr~^ 

Us = —7.16 ± 1.28 mas yr~"^ 

TT = 0.61 ± 0.22 mas, 

corresponding to a distance of 1.65^0 44 kpc. This is in agreement with the recent distances 
measured for sources related with the Cygnus-X complex (Zhang et al. 2012; Rygl et al. 
2012) and significantly larger than the value obtained by hinder et al. (2009) for Cyg 
0B2 ^5 itself. The post-fit rms in this case is 0.13 mas in right ascension and 0.55 mas 
in declination. The reduced-x^ obtained in right ascension using the errors delivered by 
JMFIT is 0.59, suggesting that no significant systematic errors remain in our data along 
that axis. In declination, however, a systematic contribution of 0.65 mas has to be added 
quadratically to the errors given by JMFIT to obtain a reduced-x^ of 1. This is not 
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Fig. 4. — Best astrometric fit to the positions obtained from the high resolution VLB A 
images. The left panel shows the right ascension as a function of time, while the right 
panel shows the declination behaviour. The filled squares show the observed positions with 
their error bars (including the systematic contributions) while the open squares indicate the 
position of the source expected from the fits at each epoch. The dotted curves are the best 
fits to the data including a parallax and uniform proper motion. 

surprising given the elongation of the source in the north-south direction. All the errors 
quoted here include this systematic contribution. 

4. Discussion 

4.1. Nature of the radio emission from the WCR(A-B) 

The emission from the WCR(A-B) has long been thought to be of non-thermal 
origin because (i) the fiat radio spectrum of Cyg OB #5 in VLA observations when this 
source is "on" compared to its positive spectral index when it is "off" suggests that the 
spectral index of the WCR(A-B) is negative, and (ii) Ortiz-Leon et al. (2011) measured a 
brightness temperature of ~ 2 x 10^ K suggestive of a non-thermal process. The brightness 
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temperatures reported here are about 10^ K, clearly favoring a non-thermal (presumably 
synchrotron) nature. 

The best-fit model of Kennedy et al. (2010) has orbital parameters (orbit inclination 
i = 90° ± 40°, argument of periastron u = 319° ± 3° and eccentricity e = 0.7 ± 0.04) that 
result in a structure as shown in Figure (the plot is made in the plane of the orbit, 
perpendicular to the plane of the sky). The large circle shows the nominal size (radius 23 
mas) of the free-free emission region from the Cyg 0B2 #5 A contact binary together with 
its uncertainty (12 mas; Rodriguez et al. 2010). It is clear that near apastron (and only 
near apastron), Cyg 0B2 ^5 B (and its associated WCR(A-B)) are located outside, or at 
least near the outer edge, of the optically-thick ionized region. This occurs because the 
orbit is quite eccentric, and explain why the non-thermal emission is detectable in spite of 
the small separation between the WCR(A-B) and the contact binary. Note that since Cyg 
0B2 7^5 B will always be located to the east of Cyg 0B2 ^5 A at apastron, the curvature 
of the WCR(A-B) is expected to always be oriented in nearly exactly the same direction, 
as observed in our VLBA data. On the other hand, we note that the decline in flux of the 
WCR(A-B) around periastron can be influenced by other plasma and radiation processes 
(Pittard & Dougherty 2006; Pittard 2009). 

As we mentioned before, the flux of the WCR(A-B) is consistent with a value of 3 
mJy from 2010 December to 2011 December. There is, however, one exception: the flux 
in 2011 January is nearly 5 times smaller than this average value. Our discussion of the 
absolute flux accuracy (end of Section 2) shows that our fluxes are accurate to about 10% 
so the decrease by a factor of 5 observed at that epoch cannot be ascribed to instrumental 
effects. Such strong variability is not expected in a wind collision region. The s=0 model 
of Kennedy et al. (2010) predicts that the flux of the WCR(A-B) is closely constant in 
the "on" state. A possible solution to this discrepancy would be for the wind driven by 
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the contact binary to be inhomogeneous. Under this assumption, two scenarios would be 
possible. The first one would be for Cyg 0B2 ^^5 B to pass through a region where the 
wind driven by Cyg 0B2 #5 A is less dense. In such a situation, the shock resulting from 
the wind collision might be weak, and insufficient to produce relativistic electrons and, 
therefore, detectable compact radio emission. Araudo et al. (2012, in preparation) computed 
the cooling time for the relativistic electrons at the WCR(A-B) considering synchrotron 
losses. Inverse Compton (IC) scattering, relativistic bremsstrahlung, and advection and 
diffusion escape of the relativistic particles from the shocked winds. They found that 
due to the dense photon stellar field provided by the unseen companion, IC scattering is 
the most important cooling channel for the relativistic electrons at energies above ~ 10® 
eV (considering different spectral types for Cyg 0B2 #5 B), implying that most of the 
relativistic electrons are cooled by IC losses. Below this energy, the cooling is dominated by 
advection of the particles on a time tesc ~ 4/bs/'i^oo ~ 4 months. Here /bs is the size of the 
WCR(A-B) and v^o the terminal velocity of the stellar winds. Since the separation between 
our successive observations is significantly shorter than this, we would not expect the fiux 
to drop appreciably between observations. A second way to produce a drop in radio fiux 
would be to have a sudden increase in the column density of ionized material along the line 
of sight to the WCR(A-B), for instance as a result of a dense clump passing in front of the 
WCR(A-B). This would result in an increase in the free- free opacity. This second scenario 
fits more naturally the observed light curve. 

4.2. The distance to Cyg OB2 #5 

As we mentioned in the Introduction, there is some significant uncertainty on the 
distance to Cyg 0B2 #5. On the one hand, recent trigonometric parallax measurements 
to masers associated with several stars in the Cygnus X star forming complex, including 
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Fig. 5. — Graphical sketch of the orbit Cyg 0B2 ^5 B (red hne). The blue circle shows the 
extent of the optically-thick photo-ionized region surrounding Cyg 0B2 #5 A (the full line 
is the nominal value, while the dashed versions show the minimum and maximum values; 
from Rodriguez et al. 2010)). As the sizes depends on the distance assumed, the values 
of both R.A. and depth are presented here in angular units (mas). The black line is the 
angular separation observed in the plane of the sky. The red squares shows the position of 
Cyg 0B2 B at three representative epochs when the WCR(A-B) can be detected to show 
the direction of the curvature of the WCR(A-B) (red arcs) in these epochs. 

the Cygnus 0B2 star forming region, suggest a distance of 1.4-1.6 kpc (Rygl et al. 2012; 
Zhang et al. 2012). On the other hand, hinder et al. (2009) obtained a distance modulus 
for Cygnus 0B2 ^^5 that yields a significantly shorter distance of 925 pc. The latter 
authors noted that theirs results were strongly dependent on the assumed stellar effective 
temperature. The value used was 36,000 K, as obtained from a spectrum analysis carried 
by Rauw et al. (1999), that suggested an 06-7 la star. They also mention that a somewhat 
unlikely effective temperature of 48000 K would be needed if the system were located to 
a distance of 1.7 kpc. The parallax obtained here, although uncertain, favors a distance 
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of order 1.6 kpc, in agreement with those obtained by Rygl et al. (2012) and Zhang et al. 
(2012) and significantly larger than the distance of 925 pc proposed by hinder et al. (2009) 

The distance to the Cyg 0B2 ^^5 system can also be estimated using an orbital 
parallax method. From our VLBA observations and the argument of periastron (u) and 
inclination (i) of the Cyg 0B2 ^5 AB orbits derived by Kennedy et al. (2010), we can 
derive the de-projected angular separation between Cyg 0B2 # 5 A and Cyg 0B2 # 5 B. 
In all four models discussed by Kennedy et al. (2010) the inclination is close to 90° (in the 
range of 85° to 90°) and can be ignored in the derivation. Then, the de-projected angular 
distance is simply equal to the projected angular distance divided by sin(360°-a;). For the 
case of the best fit model s=0, we find that at apastron the de-projected angular distance 
is 18.3 maqj. The total mass of Cyg 0B2 ^5 A was calculated to be 41.5 ± 3.4 Mq by 
hinder et al. (2009) and the mass of Cyg 0B2 ^5 B was estimated to be 23^^4 Mq by 
Kennedy et al. (2010). Thus, the appropriate mass range for the Cyg 0B2 ^5 AB system 
is between 50 and 90 Mq. From the orbital parameters given by Kennedy et al. (2010), we 
can calculate the true separation (in cm) at apastron for these two masses, and from the 
comparison between the true separation and the de-projected angular separation, deduce 
the distance to the system. We obtain distances between 1.26 kpc and 1.53 kpc. This is in 
excellent agreement with the parallax determination obtained earlier, and with the distance 
to Cyg 0B2 estimated by Rygl et al. (2012) and Zhang et al. (2012). We note, on the other 
hand, that if we assume a distance of 0.95 kpc, the total mass of the triple system would 
have to be only 21.4 Mq, lower than the primary alone. At the other extreme, assuming a 
distance of 1.7 kpc would yield a dynamical mass of 122.4 Mq, larger than usually expected. 



"^Strictly, this is the separation between Cyg 0B2 ^ 5 A and the WCR(A-B) associated 
with Cyg 0B2 7^ 5 B, but Ortiz-Leon et al. (2011) argue that Cyg 0B2 7^5B is very near 
the apex of the WCR(A-B). 
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Unfortunately, this method depends critically on the orbital parameters, and adoption of 
any other models in Kennedy et al. (2010) will give a different distance. 

Yet another (independent) distance determination can be obtained as follows. The 
separation between Cyg 0B2 ^^5 B and the WCR(A-B) can be estimated from a momentum 
ratio analysis. From Canto et al. (1996), we can relate the separation Rq between Cyg 0B2 
#5 B and the WCR(A-B) to the separation Ri between the WCR(A-B) and Cyg 0B2 #5 
A by: 

Ro = ri^'^Ri: (1) 

where t] = M^joVwo / My^iVwi is the ratio between the wind momentum ratio of the two 
stars, hinder et al. (2009) suggest a distance-independent mass loss rate of ~ 2.1 x 10~^ 
Mq year~^ for Cyg 0B2 ^5 A. Its terminal velocity was suggested to be 1500 km s~^ by 
Kennedy et al. (2010) based on the absorption of the P-Cygni Hel 1.083 fim line observed 
by hinder et al. (2009). 

The value of 1] will, of course, depend on the mass-loss rate and wind velocity of Cyg 
0B2 #5 B. In Table [3] we summarize the parameters of massive stars with luminosity class 
V. The effective temperature, luminosity, radius, spectroscopic mass, and the flux of ionizing 
photons are taken from Vacca et al. (1996). The mass loss rate was calculated following the 
prescription of Vink et al. (2001) for stars of solar metallicity, and the terminal velocity of 
the wind, Voo, was calculated assuming foo/'^esc = 2.6 (Vink et al. 2000), where fesc is the 
escape velocity from the stellar surface. Finally, we calculated the expected flux density at 
8.46 GHz assuming a fully ionized, pure hydrogen wind at a distance of 1 kpc and following 
the formulation of Panagia & Felli (1975). We also used Equation 1 of Rodriguez & Canto 
(1983) to verify that the ionizing photon rate of the stars was sufficient to fully ionize their 
respective winds. 
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Table 3: Parameters of massive stars with luminosity class V. The 4.8 GHz flux is estimated 
for a distance of 1 kpc. Additionally, the value of rj for each type of star is hsted in the 
rightmost column, assuming that its wind would interact with that of Cyg 0B2 #5 A. 
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log(M) 
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1.55 
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2800 
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2700 


0.044 
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0.97 
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0.029 
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08.5 
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0.95 


1.37 


48.72 


-6.41 
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0.017 


0.032 
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5.061 
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1.35 


48.56 


-6.61 


2600 


0.009 


0.020 
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34620 
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0.93 


1.32 


48.38 


-6.79 


2500 


0.006 


0.013 


BO 


33340 


4.881 


0.92 


1.29 


48.16 


-6.99 


2500 


0.003 


0.008 
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32060 


4.789 


0.90 


1.26 


47.90 


-7.21 


2400 


0.002 


0.004 
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Cyg 0B2 ^5 B must drive a wind strong enough to produce the WCR(A-B) when it 
interacts with the winds of the echpsing binary Cyg 0B2 #5 A. As discussed by Kennedy 
et al. (2010) and Ortiz-Leon et al (2011), this suggests that Cyg 0B2 #5 B is of BO. 5 or 
earher spectral type. The value of r] for a BO. 5 star (Table [3]) yields an upper limit on the 
distance to Cyg 0B2 #5 of 1.44 kpc. This result is, again, in agreement with our measured 
trigonometric parallax, and with distances proposed by Hanson (2003), Rygl et al. (2012) 
and Zhang et al. (2012). 

Running the same argument in reverse, we can place an upper limit on the value of t] 
considering the highest possible mass for the system (90 Mq) and shortest possible distance 
(1.26 kpc). This leads to ?7 = 0.047, corresponding to a spectral type 08 for Cyg 0B2 ^5 
B. This explains a posteriori why Cyg 0B2 #5 B contributed negligibly to the thermal 
radio emission of the system: the total free-free emission from an 08 star at the distance of 
Cyg 0B2 #5 being only about 30 /iJy. 

It is important to keep in mind that our results depend strongly on the orbital 
parameters reported by Kennedy et al. (2010). Fortunately, once a model is selected, their 
largest errors are in the inclination of the orbit, while our calculations depend most strongly 
on the argument of periastron and the eccentricity. 

5. Conclusions 

In this paper, we presented an analysis of a series of multi-epoch VLBA observations of 
the compact wind collision region in the Cygnus 0B2 ^5 quadruple system. The brightness 
temperatures for the most compact emission is ~ 10^ K, clearly suggestive of a non-thermal 
process. The total measured flux is consistent with a constant value of ~ 3 mJy for the 
first seven epochs, but falls below the detection limit for the latest observations. This is 
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consistent with the interpretation proposed by Kennedy et al. (2010) that the shocked 
region plunges into the photo-ionized region surrounding the primary in the system. In 
addition, in one of the early epochs, the flux drops signiflcantly. We suggest that for this 
epoch, a dense clump of the inhomogeneous wind intersect the sight of view, and the 
free-free emission partially hides the non-thermal emission from the wind collision region. 

The measured trigonometric parallax for the system corresponds to a distance of 
1.651q;44 kpc. In addition, the distance to the system was estimated using two indirect 
methods, both yielding distances of order 1.3-1.4 kpc. These results are in agreement with 
recently measured trigonometrical parallaxes to masers related with objects in the Cygnus 
X star forming complex, and discard the small value of 950 pc suggested by Linder et al. 
(2009). 

A new set of VLBA observations when the WCR(A-B) comes back to its "on" state, 
combined with the positions presented in this work, could be used for a better determination 
of the trigonometrical parallax of the Cygnus 0B2 ^5 quadruple system. We note that in 
agreement with Kennedy et al. (2010) the next "on" state will be on 2014 summer. 
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Stiftung for financial support. The National Radio Astronomy Observatory is a facihty 
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